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ABSTRACT 

We present optical spectra of (24) Themis and (65) Cybele, two large main- 
belt asteroids on which exposed water ice has recently been reported. No emission 
lines, expected from resonance fluorescence in gas sublimated from the ice, were 
detected. Derived limits to the production rates of water are < 400 kg s~^ (Sc), 
for each object, assuming a cometary H2O/CN ratio. We rule out models in 
which a large fraction of the surface is occupied by high albedo ("fresh") water 
ice because the measured albedos of Themis and Cybele are low (~0.05 - 0.07). 
We also rule out models in which a large fraction of the surface is occupied 
by low albedo ("dirty") water ice because dirty ice would be warm, and would 
sublimate strongly enough for gaseous products to have been detected. If ice 
exists on these bodies it must be relatively clean (albedo >0.3) and confined to 
a fraction of the Earth-facing surface <10%. By analogy with impacted asteroid 
(596) Scheila, we propose an impact excavation scenario, in which 10 m scale 
projectiles have exposed buried ice. If the ice is even more reflective (albedo > 
0.6) then the timescale for sublimation of an optically thick layer can rival the 
~10^ yr interval between impacts with bodies this size. In this sense, exposure 
by impact may be a quasi steady-state feature of ice-containing asteroids at 3 
AU. 

Subject headings: minor planets, asteroids; comets: general; comets: individual 
(Themis, Cybele); solar system: formation 



1. Introduction 



Exposed water ice on main-belt asteroids is thermodynamically unstable, on timescales 
that are very short compared to the age of the Solar system. Ice is therefore not expected 
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on the surfaces of asteroids, although it may be present underneath (Fanale and Salvail 
1989). Despite this, a growing number of observations suggest that ice does indeed exist at 
or very close to (within a thermal skin depth of) the physical surfaces of some main-belt 
asteroids. For example, repetitive activity observed in two kilometer-scale, so-called "main- 
belt comets" at heliocentric distances ~3 AU (133P/Elst-Pizarro and 238P/Read) is best 
explained by the sublimation of near surface ice (Hsieh et al. 2010, 2011). However, no gas 
has been detected in these objects and they are too small (few km) and faint to allow the 
detection of water ice in near infrared reflection spectra (Licandro et al. 2011a, Rousselot et 
al. 2011). Conversely, the 100- kilometer-scale asteroids (24) Themis and (65) Cybele show 
broad, near infra-red absorptions at ~3.1 /im that are compatible in position and shape with 
vibrational bands in crystalline water ice (Campins ct al. 2010, Rivkin ct al. 2010, Licandro 
et al. 2011b). These large, low (few percent) geometric albedo asteroids raise questions 
about the stability and supply of ice in the asteroid belt, and about the possible connections 
between large, apparently ice-covcrcd but inactive asteroids and small, mass-losing main-belt 
comets in which ice is suspected but undetected. 

An alternative explanation for the 3.1 /xm band in terms of the mineral goethite (FeO(OH)) 
has been proposed by Beck et al. (2011). One difficulty with this interpretation is that, when 
found in meteorites, goethite is a product of aqueous alteration in the terrestrial environment. 
Extraterrestrial goethite in freshly-fallen meteorites is unknown (Alan Rubin, personal com- 
munication, 2011 July 01). This could mean that goethite-rich asteroids exist but that there 
is no dynamical pathway for their coUisionally produced fragments to reach the Earth as me- 
teorites or, more simply, that there is no goethite in the asteroids. (We note that the closely 
related mineral magnetite (Fe304) is found both in meteorites and in the spectra of asteroids 
(Yang and Jewitt 2010). It is not obvious how magnetite-containing meteorites could find 
their way to Earth while goethite-containing ones could not). Regardless, the interpretation 
of the 3.1 yum band in Themis and Cybele as due to surface ice remains plausible. 

The bulk density of C-type asteroids with diameters in the 100 km - 300 km range is p 
= 2000 kg m~^, albeit with a large uncertainty (Baer and Chcslcy 2008). With this density, 
the gravitational escape speeds from Themis and Cybele arc = 105 m s^^ and 140 m s^^, 
respectively. In both cases, is small compared to the ~750 m s~^ thermal speeds of water 
molecules at temperatures (~180 K) appropriate to low albedo asteroids in the outer belt, so 
that sublimated molecules should flow unimpeded into space. This leads us to ask whether 
evidence might be found for gas produced by sublimation of the ice from these bodies. In 
this short paper, we describe spectral line observations taken in search of this gas. 
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2. Observations 

The strongest limits to the production of gas in comets are placed using optical spectra 
of daughter molecules (fragments of stable parent molecules produced by photodissociation). 
The brightest daughter emission line in the optical spectra of comets is the OH 3080A band. 
Unfortunately, this band lies near the ozone atmospheric absorption edge making ground- 
based observations of OH consequently difficult, while no suitably sensitive space-based 
instruments exist. Spectroscopic detections of gas in comets are more usually made using 
the second-brightest line, produced by resonance fluorescence of the CN radical and located 
near 3889A (Schleicher 2010). Whereas OH is produced as a photodissociation product of 
H2O, the source of CN is less certain. It is at least in part a photodissociation product of 
HCN, but may also be released from solid grains (e.g. see recent discussion in Paganini et 
al. 2010). 

Accordingly, we took spectroscopic observations targeting outgassed CN from Themis 
and Cybele, on UT 2010 September 11 at the Keck I 10-meter telescope atop Mauna Kea, 
Hawaii. We used the Low Resolution Imaging Spectrometer (LRIS, Oke et al. 1995), which 
is equipped with two 4096x4096 pixel charge-coupled device detectors producing an image 
scale 0.135" per pixel and fed by a dichroic filter. We employed the "460" dichroic, which 
has 50% transmission at 4900A, and a 1" wide slit for all observations. The geometrical 
circumstances of the observations are given in Table ([T]) while the orbital properties of 
Cybele and Themis are listed in Table (|2]). 

Themis and Cybele were identified using finder charts and confirmed by their motion 
relative to the fixed stars. Integrations of 400 s and 330 s were secured on Themis and 
Cybele, respectively, together with spectra of nearby solar analogs HD209847 and Hyades 
64. During all observations we used an image rotator to hold the long axis of the 1.0" 
wide slit perpendicular to the local horizon, thereby minimizing effects due to atmospheric 
dispersion. We guided the telescope by hand using reflection from the spectrograph slit jaws 
to keep the target objects correctly positioned within the slit. For each new pointing of 
the telescope, we secured wavelength and flat-fleld calibration spectra using a set of lamps 
internal to LRIS. 

Data reduction was performed by flattening the wavelength-calibrated spectra and ex- 
tracting the signal from a region along the slit 6.75" in length. Reflectivity spectra were 
computed by dividing each object spectrum by the spectrum of a solar analogue and nor- 
malizing the result to unity at 3900A. The normalized reflectivity spectra are shown in 
Figure ([1]), where they are offset vertically for clarity of presentation. Themis and Cybele 
are spectrally featureless at the wavelengths of CN. 
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We estimate a limit to the gas production rate using the procedure described in Jewitt 
et al. (2009). An upper hmit to the CN flux density, fcN, is obtained from 



fcN = 5a, 
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Here, cxc is the (dimensionless) uncertainty on the normalized continuum at the wavelength 
of the CN band, fs and /sggo (erg cm~^ s^^ A^^) are the flux densities at the wavelengths 
of the B-filter and at 3880 A, respectively, and Wcn =40 A is the width of the portion of 
the spectrum in which the CN emission is sought. Quantity S^nm/SB is the ratio of the 
reflectivities at the wavelengths of CN and the B-filter. Both asteroids are nearly spectrally 
neutral (Table [3]) and we set S'sgso/'S's = 1- Quantity {fssso/ fs)® = 0.52 is the ratio of the 
flux densities at the two specified wavelengths in the Solar spectrum (Arvesen et al. 1969). 
The factor of five corresponds to our deliberately conservative use of a 5a confidence limit. 

We first computed dc, the standard deviation on the reflectivity in two continuum 
windows flanking the CN band. We used the continuum windows Ci = 3800 - 3870A and 
C2 = 3910 - 4OOOA and omitted the wavelengths of the solar H and K lines of Calcium 
since these were imperfectly cancelled by the computation of the reflectivity spectrum. To 
measure Uc we first removed residual trends by fitting a 2nd order polynomial to the 3800 A < 
A < 4000 A region, excluding wavelengths from 3870 - 39 10 A and divided the reflectivity 
by the polynomial. For Themis and Cybele, we obtained o"c = 0.0026 and (Xc = 0.0027, 
respectively. 

We determined fs by a self-calibration method, as follows. The 1" slit excludes a signif- 
icant but unknown fraction of the light from the asteroid, and so the signal in the spectrum 
cannot be used directly to secure accurate photometry with reference to the standard stars. 
However, the fraction of the total light excluded by the slit is approximately independent 
of wavelength, since the point-spread-function of LRIS is constant and since the use of an 
image rotator precludes the possibility of differential refraction out of the slit. Given this, 
we use the flux density calculated from the broadband magnitudes (Table [3]) to estimate 
the flux density at the wavelength of the CN band. This method is subject to errors in the 
published B-filter magnitudes owing to, for example, possible rotational lightcurves of the 
objects. However, both objects have small rotational variations (~0.09 magnitudes with a 
rotation period of 8.3744±0.0002 hr for Themis (Harris et al. 1989) and ~0.06 magnitudes 
and 6.07±0.05 hr for Cybele (Schober et al. 1980)) and in neither case is rotation a likely 
cause of significant photometric error. As a result, we feel that the uncertainties in fs in 
this self-calibration method are probably no more than a few tens of percent, and certainly 
smaller than a factor of two, which is sufficient for our present purposes. The final 5a limits 
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to the fluxes from Tliemis and Cybele are fcN = 7.5x10 erg cm ^ s ^ and 6.1x10 erg 
cm~^ s~^, respectively. 

Tlie CN band is produced by resonance fluorescence of sunliglit, in wliicli cometary 
molecules scatter photons randomly in direction but without change in energy. When the 
coma is optically thin, the scattered flux, fcN [erg cm~^ s~^], received at distance A [cm] is 
simply proportional to the number of CN molecules, N, in the section of the slit from which 
the spectrum was extracted, or 

The constant of proportionality in Equation ([2]) is called the resonance fluorescence efficiency, 
g{R) [erg s~^ molecule"^]. It is conventionally evaluated at heliocentric distance R = 1 AU 
and scaled to other distances using g{R) = g{l)/R^. At a given R, the fluorescence efficiency 
of CN can vary by a factor of ~1.5, owing to the Doppler shift between cometary lines and 
the Solar spectrum (the so-called Swings effect). On UT 2010 Sep 11, Themis and Cybele 
had radial velocities dR/dt = +0.7 km s^^ and +1.1 km s~^, for which the fluorescence 
efficiencies are, respectively, g{l) = 2.8x10^^^ [erg s~^ molecule"^] and g{l) = 2.9x10^^'^ 
[erg s~^ molecule"^] (Schleicher 2010). 

Substituting into Equations ([1]) and ([2]), with R and A drawn from Table [2] for UT 2010 
Sep 11, we obtain = 7.7x10^^ CN molecules in the projected slit at Themis and A^ = 
8.1x10^'' CN molecules at Cybele. To estimate the production rate from emission line data 
we used the Haser model, in which the parent of CN decays with a length scale (at R = 
1 AU) ip = 1.3x10^ km and the CN daughter itself is photo-destroyed on the length scale 
id = 2.1x10^ km (A'Hearn et al. 1995). We integrated the Haser number density model 
over the projected spectrograph slit and assumed that the outflow speed of the sublimated 
gas is V = 0.75 km s~^, consistent with measurements in comet C/Hale-Bopp at about this 
heliocentric distance (Biver et al. 2002). For Themis, the resulting limit to the production 
rate of CN is Qcn < 3.6x10^^ s"^ For Cybele, the corresponding limit is Qcn < 3.5x10^^ 

We are interested in the mass-dominant volatile H2O, not CN. In other comets, the 
average ratios of the production rates of CN to OH and of OH to H2O are Q oh /Qcn ~ 320 
and Qoh/Qh20 = 0.9, respectively, giving Qh2o/Qcn ~ 360 (A'Hearn et al. 1995). If this 
ratio applies to Themis and Cybele, the implied limits to the mass production rates in water 
are QH2O < 1.3x10^^ s~\ corresponding to < 390 kg s~^ for Themis and QH2O < 1.2x10^* 
s~^, corresponding to < 370 kg s~^ for Cybele. The accuracy of these limits is subject 
to the assumption that Qh2o/Qcn has a value similar to that measured in comets. We 
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note, however, that Lovell et aL (2010) independently used UV and radio-hne observations 
to directly obtain an upper limit Qoh < 1x10^^ s~^ (3cr) at Themis, corresponding to 
300 kg s^^. This value is consistent with the limits derived here based on CN. 



We seek to estimate the fraction of the surface on each asteroid that could be occupied 
by exposed water ice, reasoning as follows. The upper limit to the mass loss rate from 
our spectra sets an upper limit to the area of exposed ice on the surface. The limit is a 
function of the unknown albedo of the ice, with dirty (low albedo, warm) ice confined to a 
smaller possible area than clean (high albedo, cold) ice. Separately, the mean albedos of the 
asteroids set an independent albedo-dependent constraint on the area of exposed ice, since 
a large area of high-albedo ice, even when averaged with the lower albedo of surrounding 
non-ice surface, would violate the low measured albedos. The two observational constraints, 
on mass loss and on the average albedo, limit both the possible ice albedo and the fraction 
of the surface covered by ice. 

If a fraction of the surface, /, is covered in ice of albedo pi and the remaining fraction, 
1 — /, is covered in non-ice (or "dirt") of albedo pd, we may write the average albedo as 
P = [fPi + (1 - f)Pd]- Solving for / gives 



Two black lines in Figure ([2]) show Equation (jS]) with two assumed values of the non-ice 
albedo, pd = 0.02 and 0.04, to represent the least reflective asteroid surfaces known. Accept- 
able solutions must lie along or between the trajectories defined by these curves in order to 
guarantee that the area- averaged albedo equals the measured value (from Equation Q). 

To compute a second constraint from the absence of detected gas, we solved the heat 
balance equation on the surface of a spherical model asteroid having the instantaneous 
location of Themis and Cybele as in Table ([1]). The equilibrium sublimation rate is given by 
solution of 



3. 



Discussion 



3.1. Simple Model 
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in which Lq is the luminosity of the Sun, R is the hehocentric distance in meters, A is the 
Bond albedo, 9 is the angle between the direction to the Sun and the normal to the surface, e 
is the emissivity of the body, a is the Stefan-Boltzmann constant and L{T) is the latent heat 
of sublimation of water ice at temperature, T. The quantity dm/dt is the sought-after mass 
flux of sublimated ice. The two terms on the right represent power radiated from the surface 
into space and power used to break hydrogen bonds in the sublimation of ice. The term on 
the left represents power absorbed from the Sun. Equation (jl]) ignores a term to account 
for conduction of heat into the interior, but this neglect is justified (in the case of Cybele, 
at least) by detailed thermophysical models in which the thermal inertia is found to be a 
very small / = 15 J m~^ s"'^'^ (Mueller and Blommaert 2004). This low / corresponds 
to a highly porous surface material in which heat conduction will be negligible. Equation 
dl]) cannot be solved explicitly, since it contains two unknowns, namely T and dm/dt. We 
used experimental determinations of the temperature dependence of L{T) and of dm/dt from 
Washburn (1926) and solved Equation numerically for each element of the surface on a 
sphere located at the distance of each asteroid. Figure ([3]) shows sample solutions for dm/dt 
as a function of angular distance from the subsolar point, on a model nucleus in which the 
spin vector points at the Sun. Again, in the case of Cybele, a sun-pointing spin vector is 
justified by a detailed analysis. The geocentric ecliptic latitude and longitude of Cybele on 
2010 Sep 11 were A = 53°, /3=-4°, respectively, which lies within 21±15° of the pole solution 
reported by Mueller and Blommaert 2004. The spin vector of Themis is unknown, but time- 
resolved spectra over one night show no temporal variation (Campins et al. 2011), consistent 
with a nearly pole-on perspective. 

Figure ([3]) shows that the subsolar, thermal equilibrium sublimation flux for A = is 
about 10~^ kg m~^ s~^, but this falls by ~4.5 orders of magnitude as the albedo rises to A 
= 0.8 and also drops precipitously at 9 >80°. The stability of the ice evidently depends very 
strongly on whether it is "clean" (high albedo) or "dirty" (low albedo) and on the position of 
the ice relative to the subsolar point (Figure [3]). The difference between dirty and clean ice 
can, in this context, be the result of the admixture of very small mass fractions of absorbing 
material. For example, <C 1% contamination by mass can, depending on the grain size and 
composition of the contaminant, depress the optical albedo of ice to <0.1 (Clark 1982). 

We computed the maximum possible ice-covered surface fraction by integrating dm/dt 
over the coldest areas on the model asteroid, from / = (1 — sin(6'c)), with 9c defined by 



In Equation ([5]), dM/dt (kg s ^) is the total mass production rate (which we set equal to the 




(5) 
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upper limit values obtained from the data) and 6c is the critical subsolar angle, such that 
sublimation integrated over all regions with 6 > 6c equals dM/dt. The asteroid radius, r„, 
is taken from Table ([3]). 

Figure ([2]) shows the sublimation constraints graphically in terms of the fraction of the 
surface occupied by ice as a function of the albedo of the ice. Red and blue curves show the 
constraints imposed by the non-detection of sublimated molecules for Themis and Cybele, as 
determined by Equations (jl]) and Acceptable solutions must lie on or below the plotted 
red and blue curves, otherwise gas production by sublimation would have been large enough 
to be detected. The region satisfying both constraints is shown shaded. Figure ([2]) shows 
that models in which ice is globally distributed (/ = 1) on Themis and Cybele are excluded. 
At one extreme, this is because full-surface clean ice would violate the (low) average albedos. 
At the other, full-surface dirty ice would be warm enough to strongly sublimate and would 
violate the spectroscopic constraint on sublimated gas. No intermediate albedo solution 
with / ~ 1 is found. Instead, the allowed solutions for both asteroids are those with ice 
albedos pi > 0.3 and coverage fractions / < 0.1. Ice on these bodies, if it is present, must 
be relatively clean and spatially localized. 

Measurements from a single rotation show that the ice band on Themis is approximately 
constant in depth (Campins et al. 2010), consistent with a small angle between the line of 
sight and the spin vector. Independent observations dispersed over a five year interval show 
a larger variation (e.g. using spectra reported in the Supplement to Rivkin and Emery 2010 
we measured the fractional band depth in 2005 to be ~40% deeper than in 2003 (continuum 
linearly interpolated between 2.7 and 3.5 /xm). Some of this variation is no doubt instrumen- 
tal, but the data are consistent with a substantial variation). In the future, a determination 
of the spin vector of Themis will be of great value in further localizing the ice on this body. 
Spin vector solutions have been published for Cybele (Muller and Blommaert 2004) but 3 
yum spectra have been reported from only two nights (UT 2009 Sep. 9; Licandro et al. 2011), 
insufficient to determine the spatial distribution. 

It is interesting to note that the 3.1 yum band appears in Themis and Cybele at only 
~10% of the local continuum depth. This band is intrinsically very strong (absorption coeffi- 
cient k ~ 10^ m~^; Irvine and Pollack 1968) and frequently appears saturated in astronomical 
spectra. The subdued band depth is consistent with a 10:1 dilution of the ice spectrum with 
bland (dust) continuum and so matches the inference drawn above, that the ice surface cov- 
erage fraction is / < 0.1. Rivkin and Emery (2010) took a different approach, arguing that 
the unsaturated band results from a very limited optical path length through the ice. They 
used ice coatings only ~45 nm in thickness to model the spectrum of Themis, apparently 
assuming / = 1, and were able to match the band depth and shape in detail. The authors 
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speculated that a thin, global ice film could result from freezing of water vapor emanating 
from beneath the surface. 

However, the thin, global ice film hypothesis has two problems. First, the radial tem- 
perature gradients on the sun-facing hemispheres of Themis and Cybele are in the wrong 
direction for surface frost to form. The dayside surfaces are warmer than the subsurface 
regions so that any free water molecules would migrate towards, and freeze at, deeper layers, 
not at the hot surface. Second, even if a 45 nm ice layer could form (by freeze out on the night 
side, for example), its lifetime to sublimation would be very short, requiring an unreasonably 
large water mass flux from each asteroid in order to maintain steady state. To see this, we 
first note that a 45 nm ice coating is too thin to have a substantial effect on the albedo, as 
shown by both models and measurements of non-absorbing coatings on grains (e.g. Lasue 
et al. 2007). As a result, the temperature of the ice film would be set by the temperature of 
the grain material with which it sits in close physical contact. For low Themis-like albedos. 
Figure ([3]) shows that the subsolar sublimation rate is of order dm/dt = 10~^ kg m^^ s^^, 
falling to dm/dt = 10~^ kg m~^ s~^ at a subsolar angle of 60°. At these rates, an ice film 
45 nm in thickness would sublimate away on timescales of only 5 s to 50 s over most of the 
sunlit surface of the asteroid. Rivkin et al. (2010) and Campins et al. (2010) assert that the 
ice is widespread or global, meaning that the total mass loss rate due to sublimation from 
an r„ = 100 km radius Themis or Cybele-like asteroid should be Airr'^dm/dt = 10^ to 10^ 
kg s~^, in steady state. These rates are two to three orders of magnitude larger than our 
spectroscopic limits to the gas production (<390 kg s~^ and <370 kg s^^ for Themis and 
Cybele, respectively) and so are observationally ruled out. Therefore, we conclude that the 
thin, global ice film hypothesis is inconsistent with the data. Another explanation is needed. 



3.2. Excavation Model 

Themis and Cybele are too small to have retained primordial heat, and too small to 
be significantly heated by the decay of long-lived radioactive elements trapped within their 
constituent rocks. Unlike some planetary satellites (for example, Enceladus) they also lack 
any source of fiexural heating driven by external periodic gravitational forces. Therefore, we 
reject the possibility that surface ice could represent water vapor frozen on the surface after 
being driven out from a warm interior by internally-driven fumarolic activity. 

Recent observations of another large asteroid suggest a different possibility. The 113 
km diameter (596) Scheila displayed transient activity after being struck by a body several 
tens of meters in diameter (Jewitt et al. 2011, Ishiguro et al. 2011 while Bodewits et al. 2011 
deduced a somewhat larger impactor). Dust ejected from Scheila was observed to disperse 
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under the action of solar radiation pressure on timescales of weeks, while the near infrared 
ice bands were not detected (Yang and Hsieh 2011). Wc suggest that Themis and Cybele 
are analogs of (596) Scheila in that they were also struck by small projectiles but differ from 
Scheila in that the escaping ejecta has long-since departed and that the impacts excavated ice. 
In all three cases, impact ejecta was then deposited on the surface as fallback debris, perhaps 
creating ray-like systems on each body. Sub-surface ice can survive against sublimation at 
~3 AU when insulated by even meter-thick layers of refractory regolith (Schorghofer 2008). 
As we will discuss momentarily, craters of the size envisioned for Themis and Cybele would 
easily puncture meter-thick refractory surface layers. It is thus reasonable to conjecture that 
some impacts might excavate buried ice in this way. 

We first need to estimate the size of the projectile needed to expose ice sufficient to 
cover ~10% of the surface. We assume a reference thickness of the impact deposited ice as 
£ ~ 1 mm, since an ice or frost layer much thinner than this will not strongly affect the 
albedo. Then, the volume of the exposed ice is F ~ /^rr^/, where r„ is the asteroid radius. 
Most of the ejecta from impact at 5 km s~^ onto a 100 km scale body is launched below 
the escape speed (Housen and Holsapple 2011). We therefore assume that the ice volume 
is comparable to the crater volume, Vc = 27rr^/5, where we approximate the crater by a 
cylinder with radius Tc and depth to diameter ratio 1/5. Setting V = Vc, we find 




for the radius of the impact crater needed to produce enough ejecta to cover a fraction, /, 
of the surface with ice to a depth, i. Substituting / = 0.1, r„ ~ 100 km and expressing i in 
mm, wc find r,. ~ 130 [£/(lmm)]^/^ m. For i = 1 mm, such a crater would occupy ^10~^ of 
the surface area of a 100 km radius target body. 

In the gravity-scaling limit, the radius of the projectile needed to form a crater with Tc 
= 130 m on a 100 km radius asteroid at impact speed 5 km s~^ is ~1 m, assuming normal 
impact and identical target and projectile densities. However, experiments with impact into 
porous, fragmental and icy materials show that a larger projectile is needed to form a crater 
of a given size, relative to the gravity-scaled estimate (e.g. Holsapple 1993, Arakawa et 
al. 2000). For example. Equation (1) of the latter paper gives a projectile radius ~5 m for 
normal impacts into porous, icy material, with other parameters held the same. The surface 
conditions on Themis and Cybele are unknown but it is quite likely that these bodies possess 
a fragmental surface layer caused by past impacts and statistically likely that the impacts 
were not perpendicular to the local surface, so that the projectiles were larger than indicated 
by the gravity-scaling estimate. 
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To proceed, we consider a nominal impactor radius of 10 m, on the understanding that 
this size is probably slightly too large but is in any case uncertain by a factor of several 
as a result of the unknown impactor speed, density, incident angle and impact parameter. 
We used the asteroid collision probabilities from Bottke et al. (1994) in order to estimate 
Tc, the interval between impacts of 10 m projectiles onto a 100 km radius target asteroid. 
These collision probabilities are based on measurements of asteroids larger than 1 km in size, 
requiring a factor of 100 extrapolation to reach the 10 m scale of the projectiles implicated 
here. Unfortunately, knowledge of the size distribution of sub-kilometer asteroids in the 
main-belt is limited because most such bodies are too faint to be directly observed. We rely 
on crater counts from asteroid Gaspra to provide an indication of the size distribution of 
small projectiles. There, impact craters from 0.4 km to 1.5 km in diameter are distributed as 
a power-law with a differential size index -3.7±0.5 (Belton et al. 1992, while fresh craters may 
follow a slightly steeper index according to Chapman et al. 1996). Craters in this size range 
result from the impact of projectiles a few decameters in radius, again subject to significant 
uncertainties about the surface physical properties of Gaspra. Using this size distribution, we 
estimate that the timescale for impact onto Themis or Cybele-sized asteroids is Tc ~10^ yr. 
This timescale is clearly very uncertain because of the large extrapolation from the relatively 
well-sampled asteroids at kilometer size-scales down to the 10 m sizes of the projectiles. For 
example, the estimated numbers of asteroids with r„ > 5 m range from ~10^'^ to ~10^^ 
(Davis et al. 2002), indicating a two order-of-magnitude uncertainty in at small projectile 
sizes. With these caveats in mind, we consider the implications of a 10^ yr timescale for 
impact. 

First, we note that the number of non-overlapping craters that can be placed on the 
surface is roughly Nc ~ 4r^/(r^). With r„ = 100 km and Vc = 130 m, we find Nc ~ 2x10^, 
corresponding to the accumulation of impacts over ~2 Gyr. Buried ice could persist against 
repeated impact excavation over a large fraction of the age of the solar system. 

We next determine the conditions which must prevail for excavated ice to survive on the 
surface for Tc ~ 10^ yr. The mass column density of a 1 mm thick ice layer is pi = 1 kg m~^. 
If this layer is to survive for 10^ yr (3x10^° s), the mean loss rate cannot exceed pi/Tc ~ 
3xl0~^^ kg m~^ s~^, setting an upper limit to the equilibrium sublimation temperature 
Tc < 131 K from Equation (j3]). In the specific model considered here, with the Sun on the 
projected rotation axis, sufficiently low temperatures are found in a thin band around the 
equator where cos (6') in Equation (jl]) is suitably small. By analogy with the Moon, surface 
ice on Themis and Cybele might survive best in regions protected from Sunlight by local 
topography. 
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3.3. Numerical Thermal Model 

To explore the stability of surface ice in more detail, we used a numerical (three- 
dimensional) heat transport model (Guilbert-Lepoutre et al. 2011) and considered a rotating 
body having finite thermal conductivity, a realistic rotation period (6 hr) and three obliq- 
uities 1^ = 0°, 45° and 90°. The surface thermal inertia was taken to be / = 15 MKS units 
(Mueller and Blommaert 2004), emissivity e = 0.9 and we assumed that the spin vector at 
perihelion occupied a plane that also contained the Sun. Ice Bond albedos from 0.2 < A < 
0.8 were considered. We computed the temperature on the surface as a function of position 
and time. Sample solutions for the peak annual temperature as a function of latitude are 
shown in Figure (jl]). The model results plotted there refer to a spherical body with 90° 
obliquity moving in the orbit of Cybele and with four values of the ice albedo, as marked. 
For this large obliquity, the equator shows the lowest temperatures. For albedos >0.6, the 
equatorial regions never rise above 131 K, the critical temperature for retaining surface ice 
over the 10'^ yr collision interval. 

In Figure ([5]) we show the cumulative fraction of the surface whose peak temperature 
remains less than a given value, T. For simplicity, we show only the results for albedo = 
0.6 and tabulate results for other albedos in Table (jl]). A vertical dashed line in Figure 
(jS]) shows the critical temperature T = 131 K (c.f. Figure (jl])), below which surface ice 
can survive against sublimation for longer than the 10^ yr collision time. For the particular 
case shown, ai ip = 0°, almost 50% of the surface remains at T < 131 K all around the 
orbit, a result of permanently cold polar regions illuminated obliquely. Table (jlj) shows that, 
even for A = 0.2, these polar caps survive over a substantial fraction of the surface (/ ~ 
0.28). At high obliquity, ip = 90°, the cold region shrinks to ~20% of the surface and is 
distributed instead as an equatorial band. The high obliquity models show the highest of 
all temperatures (at the poles) but retain cold nearly equatorial bands, as shown. For A > 
0.6, these equatorial bands remain below 131 K throughout the orbit and so are potential 
regions for the preservation of surface ice on 10^ yr timescales. The smallest temperature 
variation over the surface is found for intermediate obliquities {ip = 45°). For A < 0.6 and 
this obliquity, no part of the surface is cold enough to retain ice for 10^ yr (Figure [5]). 

Taken together, the numerical model results show that ice can persist for ~10'^ yr 
timescales on the surfaces of asteroids in Cybele-like and Themis-like orbits, over a wide 
range of ice albedos and asteroid obliquities. Small obliquities strongly favor the long-term 
survival of ice in cold, polar regions, but ice can also survive near the equators of asteroids 
having obliquities ~90°. For albedos A >0.6, a millimeter-thick ice layer can survive for ~10'^ 
yr on 10% or more of the surface regardless of the obliquity. In this sense, the conduction 
models are compatible with the scenario proposed above, in which buried ice is excavated 
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by one or more small impacts and spread across a fraction of the asteroid surface as fallback 
debris. While we possess no proof that this is the origin of the near infrared spectral features 
on either Themis or Cybele, the inevitability of impact and the simplicity of the impact 
excavation model provide ample motivation for an expanded search for ice-like absorptions 
in, and outgassing from, the asteroids. 
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4. Summciry 

The near infrared spectral signature of water ice has been recently reported on two 
main-belt asteroids, (24) Themis and (65) Cybele. We took optical spectra to search for 
gaseous emission from these objects caused by the sublimation of surface ice. 

1. The CN 3889 A emission line of CN was not detected in either object. Assuming comet 
-like H2O/CN mixing ratios, we compute upper limits to the water production rates 
<390 kg s~^ for Themis and <370 kg s~^ for Cybele (both are 5a hmits). These rates 
are two to three orders of magnitude smaller than outgassing rates expected from the 
thin, global ice film hypothesis of Rivkin and Emery (2010). 

2. If surface ice exists on Themis and Cybele, it must be clean (albedo >0.3) and of 
limited spatial extent (surface fraction ^10%) in order to simultaneously fit the low 
average albedos of these objects and the non-detections of gas from sublimated ice. 

3. We speculate that ice might have been exposed on the surfaces of Themis and Cybele 
by recent impacts. A projectile only 10 m in radius could blanket ~10% of the surface 
in ice to a depth of 1 mm, sufficient to produce observable infrared spectral signatures. 
Impacts onto Themis and Cybele by projectiles of this size occur on timescales ~10^ 

4. Ice on Themis and Cybele with albedo >0.6 would have a sublimation lifetime in excess 
of the 10^ yr collision time. Surface ice, if it is largely free of absorbing contaminants, 
could be a semi-permanent feature of these and other large, volatile-rich main-belt 
asteroids. 
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Table 1. Geometrical Circumstances 



Name 


UT Date 




dR/d£^ 




a 


24 Themis 
65 Cybele 


2010 Sep 11 
2010 Sep 11 


3.510 
3.696 


+0.7 
+1.1 


2.868 
3.147 


14.0 
14.2 



^Heliocentric distance [AU] 

^Ratc of change of Heliocentric distance [km s~^] 
■^Geocentric distance [AU] 
^Phase angle [°] 
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Table 2. Orbital Properties 



Name 












24 Themis 
65 Cybele 


3.129 
3.440 


0.134 
0.106 


0.76 
3.54 


2.708 
3.076 


3.548 
3.805 



^Semimajor axis [AU] 
^Orbital eccentricity 
^Orbital inclination 
^Perihelion distance [AU] 
^Aphehon distance [AU] 
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Table 3. Physical Properties 



Name 






B-V^ 


pe 


QcN^ 


dM/dt^ 


24 Themis 
65 Cybele 


198±20i' 0.067±0.012^ 
273±12i 0.050±0.005' 


13.49 
13.64 


0.64 

0.68-0.73 


8.3744±0.0002 
6.07±0.05 


<3.6xl025 
<3.5xl025 


<390 
<370 



'^Diameter in km 
"^Visual geometric albedo 
'^Apparent B-filter magnitude 
^B-V color 

^Rotation period in hr 

^imit to the CN production rate in number of molecules per second 
^imit to the H2O production rate in kg 

'^For Themis, we use the value quoted on the JPL Small-Body web site but note that Chernova 
et al. (1994) independently reported albedo 0.074 

Tor Cybele, we use measurements by Mueller and Blommaert (2004). 
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Table 4. Conduction Models^ 



A = 0.2 0.4 0.6 0.8 

0.28 0.35 0.46 1.00 

45 0.58 

90 0.21 0.49 



^Fraction of the surface having 
peak temperature T < 131 K for 
obhquities 0°, 45° and 90° and four 
values of the ice albedo. 
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Fig. 1. — Normalized reflectivity spectra of Themis and Cybele showing the two wavelength 
regions, Ci and C2, used to define the CN continuum. The spectrum of Themis has been 
displaced vertically by 0.5 for clarity. 
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Fig. 2. — Allowable values of the ice albedo are plotted against the fraction of the surface 
occupied by ice. Black lines reflect the albedo constraint (Equation [3]). Red and blue lines 
show the sublimation constraint for Cybele and Themis, respectively. The shaded region 
shows allowable solutions for the two asteroids. 
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Fig. 3. — Thermal equilibrium sublimation mass flux versus angle from the sub-solar point, 
computed for flve different values of the ice albedo, as marked. 
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Fig. 4. — Peak temperature as a function of latitude for Bond albedos A = 0.2, 0.4, 0.6, and 
0.8 on a spherical model body moving in the orbit of Cybele with obliquity ip = 90°. The 
critical peak temperature, below which exposed ice can survive for more than 10^ yr, is also 
marked. 
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Fig. 5. — Cumulative distribution of the peak temperature calculated over the surface of a 
rotating sphere following the orbit of Cybele, for three obliquities and ice albedo A = 0.6. 
The critical temperature below which exposed water ice can remain for >10'^ yr is marked 
by a vertical dashed line. 



